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ABSTRACT 

The charged particles considered in this paper are scattered by 
random fields while they propagate along the diverging lines of force 
of a spatially inhomogeneous guiding field. Their longitudinal transport 
Is described in terns of the eigenfunctions of a Sturm-Liouville operator 
which incorporates the effect of adiabatic focusing along with that of 
scattering. The relaxation times and characteristic velocities, which 
appear in this matrix formulation of the transport problem, are graphed and 
tabulated. The particle density that results from a localized impulsive 
injection is evaluated as a function of space and time for two different 
regimes. In the first regime, where focusing is relatively weak, a 
diffusive mode of propagation is dominant, but coherent modes are also 
present, and they become prominent as the intensity of focusing increases. 

In the second regime, where focusing is strong, diffusion does not occur, 
and the propagation is purely coherent. This supercoherent mode corresponds 
exactly to the so-called scatter-free propagation of kilovolt solar flare 
electrons. Moreover, diffusive propagation in the first regime offers an 
explanation of several poorly understood aspects of solar cosmic-ray events. 
On a larger scale, focused transport provides an interpretation of many 
observed characteristics of extragalactic radio sources. 

Subject headings: cosmic rays: general - hydromagnetics 
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I. INTRODUCTION 

The diffusion equation underlies much existing work on solar and 
galactic cosmic- rays. However, if particle propagation takes plsee along a 
diverging guiding field under the influence of adiabatic focusing, the 
diffusive idealization is valid only if the mean free path for scattering 
in random fields is small compared to both the scale length for spatial 
variations of the density and the scale length for spatial variations of 
the guiding field. This paper describes the deviations from diffusive 
behavior that occur when neither of these conditions is satisfied. This 
formulation of transport theory is maLhematically similar to the approach 
that I took in three previous papers (Earl 1973, Paper I, 1974a, Paper II 
and 1974b, Paper III). Consequently, these papers and the equations therein 
will be designated below by their roman numerals. Most of these references 
are to Paper II which analyzed a coherent mode which is qualitatively similar 
to the supercoherent mode that occurs when adiabatic focusing is sufficiently 
intense. The systematic effect of focusing considered here is completely 
different from the stochastic effect considered by Golds teln,Kllmas and 
Sandri (1975), which arises from small scale divergences of the random field. 

In §11, the transport problem is formulated in terms of eigenfunctions 
of an operator which incorporates both scattering by random fields and focusing 
by a spatially inhomogeneous guiding field. The behavior of these eigen- 
functions is described in §III. Focused diffusion, a transitional mode 
which spans the gap between the purely diffusive transport that occurs when 
focusing is weak and the purely coherent transport that occurs when focusing 
la strong, is discussed in §IV. The supercoherent mode is introduced in 5V, 
Thus, the theory developed in these sections not only identifies and describes 
a novel mode of particle propagation but also establishes its relationship to 


diffusion. 
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The existence of pronounced coherent effects opens up many possibilities 
for the interpretation of astrophysical phenomena. These possibilities are 
explored in an interplanetary context in §VI where the so-called scatter- 
free propagation of solar flare electrons is explained and where several 
poorly understood aspects of solar cosmic-ray events are Interpreted. In 
IVII, the structure of extragalactlc radio sources is explained in terms of 
focused transport. Here, the transport phenomena introduced in this paper 
give rise not only to the twin lobes of radio emission, which are a basic 
feature of these sources, but also to many other details of their morphology. 
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II. MATRIX FORMULATION OF TRANSPORT THEORY 

The particle distribution function f{y, z, t} evolves according to 
the equation 


3t 3z 2 3y 3y 




( 1 ) 


In vhlch the effect of adiabatic focusing is represented by the second 
term on the right hand side (Roelof 1969), which involves the scale 
length L for spatial variations of the guiding field. 


1 

L 


1 ^ 
B 3z 


( 2 ) 


In equation (1), z is distance parallel to the mean field, y is 
the cosine of the pitch angle, V is particle velocity, and t is time. 
The Fokker-Planck coefficient for pitch angle scattering will be des- 
cribed by 


(|»{y} * A|y|** ^ (1 - y2) , (3) 

in which q is the spectral index of the power law that gives the mean 
square amplitude of field fluctuations at wave number k within an 
interval dk, terms of the spectral density at a 

reference wave number k^. Here, the parameter A, 


A » 2if ^ Q (k-r, ) ^ , 

jj^2 XX 0 L * 


(4) 


can be expressed in terms of the particle rigidity R, velocity V, and 
Larmor radius r, and the spectral parameters Q , q, and k_ (Jokipii 
1966). Although the validity of the quaslllnear approach that underlies 
equation (4) has been questioned, this simple relationship is Invoked 
here for purposes of Illustration with the understanding that the 



5 


numerical values of certain parameters, which are expressed below not 
only as functions of q and A but also in general terms, may have to be 
revised when a consensus is reached on the correct treatment of pitch 
angle scattering. Such a revision would not affect the qualitative 
validity of the conclusions reached here. 

The assumption that L is constant, which will be adopted through- 
out this paper, greatly simplifies the analysis that follows, but as 
long as L does not change much within one scattering length, it does 
not significantly limit the applicability of the results. This assump- 
tion implies that the guiding field decreases exponentially, which means 
that its lines of force diverge from one another as z increases. Because 
of this divergence, the lateral area over which particles are spread 
Increases with z, for particle transport perpendicular to the field lines 
proceeds relatively slowly. Consequently, the normalization that 


corresponds to a fixed total number of particles is 

+1 4 <» , /■+“ 

Nq => “ y* J f{y, z, t} dg dz =y < 


Fq{z. t} dz (5) 


in which the exponential factor takes into account this variation in 
the area over which the isotropic density Fq Is spread. 

When equation (1) is integrated over u from -1 to +1, the scattering 
term contributes nothing because vanishes at both limits , while the 
focusing term can be integrated by parts to yield 


!!o_ ri- + iis 

at L L J 


where S is the streaming flux defined by 


S - (V/2) 


gf dg . 
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Because the expression in square brackets Is the divergence operator 
(Roelof 1969) , equation (6) expresses an important and familiar equality 
between the temporal rate of change of the density Fq and the negative 


divergence of the flux. 

Papers I, II, and III invoked expansions of f in terms of eigen- 
functions of the scattering operator which appears as the first term on 
the right hand side of equation (!) . To describe the effects of focusing, 
this paper Invokes eigenfunctions of a very similar operator which also 
includes the focusing term. Thus, the focusing eigenfunctions 
which are defined in the same spirit as the scattering eigenfunctions 


satisfy the following equation: 






d j. V Q - h2) 1:1 4. j- q = n 


( 8 ) 


where the eigenvalue (2/a^) , which replaces the scattering eigenvalue 

(2/t ) , describes the temporal decay of an anisotropy proportional to 
K 

Q . Equation (8) can also be written in the form 
1C 


d -G 
dp ® 


dQw 2 -G 

dp Oj. N 


(9) 


where the function G{p} which appears in the exponential. 


G{p} = ? / 


1-v^ 


dv = 


V p 


2-q 


( 10 ) 


L •'o AL (2-q) ’ 

is the same as the one defined by equation (111-32) except for a scale 
factor (1/L). The important parameter (V/AI.), which is the ratio of the 
scattering length (V/A) to the scale L of guiding field variations, 


characterizes the intensity of focusing. The following boundary condi- 
tion, which is to he imposed at p « +1 and at p = -1, completes the 
specification of the eigenfunctions and ensures that they are well behaved 
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♦ 


5 

dli 


K 


0 . 


( 11 ) 


In the limit L contrast to the situation 

discussed previously, where even nutubered were even functions of u 
and odd numbered were odd functions, the have no special symmetry 
when L Is finite. In Its usual sense as a description of spatial 
symmetry, the word parity is obviously inappropriate here. Nevertheless, 
this word will be used below in its mathematical sense to designate 
whether the indices of eigenfunctions are odd or even. 

Because equation (9) has the form prescribed by Sturm~Liouville 
theory, the functions form an orthogonal set in terms of which the 
distribution function can be expressed as a series expansion 

f{p, 2i t} = ^ ‘Sc 

where the factor d^ defined by 

k “ ‘/r 


converts into a normalized eigenfunction and where the coefficients 


fj^ are given by 




Qj^ f dli 


(14) 


In equations (13) and (14) , which typify the integrals that occur when 
orthogonality is invoked, exp{-G} is a weighting function which emphasizes 
the contribution from the region p < 0 where the odd function G is 
negative. This asymmetry, which becomes very pronounced when (V/AL) is 
large, makes focused transport qualitatively different from rectilinear 
transport . 
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When equation (12) is substituted in equation (1) , orthogonality 
implies that the coefficients fj, are described by a set of 
differential equations the first four of which are 


!!c 

3 t 

3 f 


+ U 


!!c 

00 3z 
3f. 


i . + _i + U i . 

at Oj^ 11 3z 

3f2 f2 3^2 

af f3 3f3 

^ ^ ^ «33 ^ 



af . 


3 f , 



-"01 

1 . 

3z 

“"02 

2 

3z 

-"03 

3 

3z 



■^12 

'^2 

“^3 

®^3 

-%1 

3z 

3z 

3z 



■"12 


-'^23 

''3 

-^02 

3z 

3z 

3z 


'^0 

■^13 

'^1 

-«23 

'^2 

-^03 

3z 

3z 

3z 


(15) 

(16) 

(17) 

(18) 


where the characteristic velocities, 

"JK - "kJ ■ Oj 

play the role of matrix elements that couple the temporal evolution of 
the coefficients to their gradients. This matrix formulation of trans- 
port theory is analogous to that derived previously (see eq. [II-IO] 
[11-13]), but it differs in that each coefficient is coupled to its own 
gradient through the diagonal velocity elements which appear on the 
left hand sides of equations (15) - (18). Moreover, the source terms 
on the right hand sides of these equations embody coupling between 
coefficients of similar parity that did not appear in previous papers 
where the characteristic velocities corresponding to and Uj ^3 were 
zero. However, in the limit L ^ the matrix elements that couple 
coefficients of opposite parity, Uq^, 11 ^ 3 , ^23’ 

same characteristic velocities, Vq 3 » V 32 t ^^23' 


Paper II. 
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The fiist eigenfunction, “ constant, is an isotropic component 
which satisfies equation (8) provided that the relaxation tine is 
Infinite. Because all of the higher order eigenfunctions are orthogonal 
to Q^, they must satisfy 



Qj, dp = 0 


( 20 ) 


but this condition does not imply that the Isotropic density associated 
with the higher-order eigenfunctions is zero, for this weighted integral 
is not the same as the unweighted average that gives the density. More 
specifically, the density is given by 


where 



f du -E <Qr>^K 


=■ hd^ 



dp 


( 21 ) 

( 22 ) 


is the average density associated with Qj^. Similarly, according to 
equation (7), the flux is 


= - E \ ‘k 


where the velocity V defined by 

1 \ 


\ 


•'-1 


p dp 


(24) 


characterizes the flux associated with Qj^. Because the expressions for 
density and flux involve all components, the situation is different from 
that in unfocused transport where the density is identical to the 
isotropic component and where, consequently, the lowest-order matrix 
equation is also the flux equation. Instead, the derivative (BF^/Bt) 
that appears in the flux equation (eq. [6]) ™st be calculated by 
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summing the quantities <Qj^> (3fj^/3t) predicted by the matrix equations 
(eqs. [15] - [18]). In the resulting expression, 


!!o 

3t 


- '<'’o>"oi *<'*!> ”u 

- t<«o>0„2 +<Qi>Bi2 


+ <Qj>noj 

+ <02>"l2 


+ 


+ 


+ 


...I 


!!o 

iz 


‘ 32 ■ o 

I 

3z 


the sums within square brackets reduce to 


(25) 


which is an identity that follows from equation (24) when the integrand 
uQjj is expanded in an eigenfunction series with the aid of equations (14) 
and (19). A second identity, which relates V to the relaxa- 

Iv 

tion time 

«1k> \ 

~^‘r • < 27 ) 

follows from 



In which the first equality results from an integration by parts, the 
second from the following substitutions (see eqs. [9] and [10]): 
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dG . .r _ V 1 “ 
du ^ L ’ 


and the third from a further integration by parts. When these identities 
are invoked in equation (25), it reduces to the flux equation, for the 
gradient terras sum to (8S/3z) while the ( ^ ^ terms sum to (S/L) • 

If some of the matrix equations are left out of this summation, the flux 
equation is not satisfied. Consequently, when the integral specified by 
equation (5) is performed upon the solutions of a truncated set of 
matrix equations, the resulting total number of particles changes with 
time. Although this deviation from proper normalization is disconcerting, 
it has minor significance as long as the real temporal evolution is 
rapid compared to the relatively slow decay artificially introduced by 


truncation. 
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III. FOCUSING EIGENFUNCTIONS AND EIGENVALUES 

Focusing eigenfunctions are analogous to the familiar orthogonal 
functions that occur in classical and quantum physics. In particular, 
when q • 1 and (V/AL) <■ 0, they reduce to Legendre polynomials. But 
it is only in such exceptional cases that analytical methods yield closed 
expressions for eigenvalues and characteristic velocities. In general, 
these parameters must be evaluated numerically. There is n:s need to 
present here the lengthy details of this evaluation, for the final results 
are sufficient to specify completely the matrix formulation of transport 
theory. Nevertheless, it is appropriate to outline briefly the method 
that gave these results. It was basically the iterative method of 
Stodola and Vianello (Hildebrand 1949, Chap. 5) supplemented by the 
procedures described by Boul.idls and Ruggiero (1944) for determining 
higher order eigenfunctions. For each value of q, scattering eigenfunctions 
for (V/AL) = 0 were calculated by iteration starting with the approximate 
eigenfunctions derived in Paper 11-51 1 =■ ' le initial trial functions. 

Then (V/AL) was incremented in small At each step, new eigen- 

functions were calculated using the eigenfunctions from the previous step 
as initial trial functions. The main objective of this section is to 
present graphs that show how these eigenfunctions and the parameters 
derived from them as specified in §II depend upon q and A. 

For weak focusing, a workable alternative to numerical methods is 
to treat the focusing term as a perturbation whose effect can be 
approximated with the aid of standard quantum mechanical perturbation 
formulae. In a preliminary version of this paper, perturbations of the 
approximate scattering eigenfunctions given in Paper II were evaluated 
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with the aid of second-order theory. Although most of the results In 
§1V were first obtained through this procedure, it was eventually 
abandoned, because an intractable divergence of the perturbation expansions 
made it Impossible to analyze the strong focusing limit discussed in §V. 
But, in the course of this exercise, several useful relationships were 
discovered. Although the perturbation approach will not be pursued, it 
does confirm the fundamental validity of these identities, which appea' 
here as numerical coincidences. 

In the discussion that follows, three specific values of the spectral 
index are given special emphasis. The first, q = 1.0, corresponds to 
the isotropic scattering considered in classical transport theory. Thus, 
the results obtained for this index illustrate the effect of focusing upon 
classical rectilinear transport. The second index, q = 1.5, approximates 
that observed for magnetic fluctuations in space. Thus, the predictions 
obtained for this index apply to the interplanetary propagation of solar 
and galactic cosmic- rays. The third value, q = 1.9, corresponds to very 
anisotropic scattering. In this situation, the eigenfunctions can be 
evaluated numerically, but they are qualitatively similar to those for 
q > 2, where the evaluation is complicated by the divergent behavior of 
the function G defined by equation (10). Thus, the results obtained for 
this index indicate how focusing affects the coherent regime discussed 
in Paper II. 

The effect of focusing, for q = 1.5, is illustrated in Figure 1 
where the focusing eigenfunctions for (V/AL) = 1 are compared with 
the scattering eigenfunctions. I'e latter functions, at the left, exhibit 
two qualities which also characterize Legendre polynomials. (See Abramo- 
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wltz and Stegun, 1964, Fig. 22.8.) First, the number of zeros between 
W ■ -1 and p “ +1 is equal to the index. Second, and R^, whose indices 
are odd numbers, are odd functions of p, while R^ and R2, whose indices 
are even, are even functions. Because the first of these features is a 
consequence of general theorems ( Courant and Hilbert, 1953, Ch. 6.), it also 
appears in the focusing eigenfunctions, at the right, where the zeros are shifted 
slightly toward smaller values of p but where their number is unchanged. 

The second symmetry, which is a characteristic of rectilinear transport, 
is less fundamental. Thus, the focusing eigenfunctions display a 
prominent asymmetry such that their absolute magnitude is generally 
larger when ;i > 0 than it is when p < 0. This asymmetry becomes very 
pronounced for large values of (V/AL) where all eigenfunctions are small 
except in the vicinity of p = +1 where they are large and positive. 

Within the scope of this paper, a comparison of observational details 
with theory is not possible. Nevertheless, it is worth noting here the 
striking similarity of the anisotropy in Figure 1 to the angular 

distributions reported for solar protons and electrons by Nielsen, 

Pomerantz and West (1975). 

In Figure 2, eigenvalue spectra are plotted as functions of (V/AL) 
for the three spectral indices mentioned above. In all cases, the 
deviation of the eigenvalue (2/AOj^) from its unperturbed value (2 /ATj^) 
increases quadratically with small values of the parameter (V/AL) and 
linearly with large values. For q = 1.5 and q = 1.9, this linear 
increase of (2/AOj^) and (Z/Ao^) is such that these eigenvalues differ by 
a small and nearly constant separation, but for q = 1.0, the increase is 
such that they become equal at (V/AL) = 11.5. This degeneracy disappears 
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ac q » 1.2. On the other hand, when focusing is absent, the spectrum for q > 1.9 
foreshadows the degenerate behavior that appears when q > 2, for (2/AOj^) is nearly 
equal to zero and is nearly equal to Because this convergence does 

not occur when (V/AL) > 0, it can be inferred that focusing removes the 
degeneracies, discussed in Paper II, that characterize the coherent 
regime. This Inference was confirmed at q = 2.0 by a detailed analysis 
which also showed that the dependence of (2/AOj^) upon (V/AL) is purely 
linear. A tendency toward this disappearance of the quadratic regime 
is evident in Figure 2. Similarly, in the strong focusing limit, the 
existence of a quadratic regime affects the intercept of the linear 
relationship but not the slope. Thus, in this limit where (V/AL) ->• «, 
the eigenvalues vary as (2/a) o (V/L) . They do not depend sensitively 
upon the parameters q and A which describe scattering. 

Figure 3 shows how the four velocities that have finite values in 
the absence of focusing, ^23’ upon (V/AL). In 

all cases, they approach zero in the strong focusing limit. For 
sufficiently anisotropic scattering, exemplified by the curves for q = 1.5 
and q = 1.9 at the right, this approach takes the form of a monotonlc 
decrease with (V/AL) which sets in at smaller values of this parameter 
and becomes more precipitous as q Increases. For Isotropic scattering, 
illustrated by the curve for q = 1.0 at the left, the velocities and 
U 22 decrease monotonically, but and 8° through maxima and minima 
before decreasing. 

The six characteristic velocities that vanish when focusing is 
absent are shown in Figure 4. For weak isotropic scattering, they 
exhibit the intricate behavior shown at the left by the curve for q = 1, 
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but as (V/AL) Uqq, and decrease slowly toward large negative 

values. Similarly, Uq 2 increase from small negative values toward 

zero, while U22 becomes large and positive. When the scattering becomes 
anisotropic, this complicated and unsymmetrical pattern simplifies drama- 
tically as Is Illustrated at the right by the curves for q = 1.5 and 1.9. 
Here, the velocities that couple coefficients of odd parity to themselves, 

and Ujj, and to each other, are positive and relatively large, 

while the velocities that couple those of even parity to themselves, 
and U22. and to each other, Uq 2. negative and relatively large. This 
pattern gives rise to the supercoherent modes discussed in §V. In the 
strong- foe using limit, the coupling between coefficients of opposite parity, 
which causes diffusive effects, becomes weak for anisotropic scattering. 

This absence of coupling has a quantitative reality that cannot be 
adequately described by the graphical representation in figures 2 and 3. 

For example, when q = 1.5 and (V/AL) = 6, U22» which is the largest of the 
velocities that couple unlike parities, is only 6% of 1)^2 » which is the 
smallest of chose that couple like parities. When (V/AL) = 10, this ratio 
decreases to 0.1%. 

Table 1 gives numerical values of the three eigenvalues and ten 
characteristic velocities. From these entries, all of the parameters 


defined below can be calculated. 
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IV. FOCUSED DIFFUSION 


In the weak focusing regime, which coincides approximately with the 
quadratic regime of Figure 2, all of the characteristic velocities are 
large enough to be significant, but the higher-order relaxation times are 
much smaller than Under these circumstances, as was discussed in 

Paper I, the coefficients and f^, which are approximately proportional 
to 0^ and o^, play a relatively minor role. Consequently, focused 
diffusion, which is the fundamental mode that occurs here, can 
be discussed in terms of a truncated set of matrix equations in 

which the two lowest order coefficients, and f^^, are retained in 
equations (15) and (16), but the small coefficients of higher order 
eigenfunctions are neglected. These equations, which are 


3f 

1 

3f, 


0 


'^0 ^^1 

u — - =* - u — ^ 

00 3z 01 3z 


“1 ^1 ^^0 
— i i + u — ±. = - u — y. 

3t 11 3z 01 3z 


(29) 


(30) 


involve three velocities, U^^, and Uj^j^, and one relaxation time 

In the discussion that follows, they will be solved with the aid of the 
methods Invoked in Paper II. However, it is worth considering first the 
result of eliminating fj^ from equations (29) and (30) 




’^*0 ”oo “c 

®00 “u* 2 ■ *“00 “11^ 5.5t ■ 0 5. 

dZ 1 


(31) 

in which the left hand side is the telegrapher's equation (eq. [11-35]). 

The terms on the right hand side, which embody the effect of focusing, 

vanish when L -► ®. Thus, in this limit, the solutions given below reduce to 
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the well known results discussed in Paper II-5III. 

To obtain these solutions, it is appropriate to express both and 
as Fourier integrals over wavenumber k 


fQ{z,t} => y dK a{<} e 


Ikz 4- lut 




/ 


d< B{(c) e 


Ikz -f lu)t 


(32) 


in which, by virtue of equations (29) and (30), the amplitudes a and B 
must satisfy 

(o) + U„ic)a + U„, tc B = 0 , (33) 

00 01 

+ (id + ~ 8 = 0 • (34) 


Because these are linear homogeneous equations, the two frequencies for 
which they have a solution must satisfy a quadratic equation, 


(u) + Uqq<)(uj + - i/Oj^) - = 0 , (35) 

which states that the determinant of the coefficients vanishes. These 

frequencies, and , which can be written in the form 

0)^ = (1/a^) + V^tc ± - (< - (36) 

where 


= ^°i 

V. - , 


(37) 

(38) 


V 

c 


K 


1 


K 


2 


- *^<"00 “ u ' ' 


u 


01 


2^1 V/ 
^1 “ «00 




( 39 ) 

( 40 ) 


( 41 ) 
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correspond to normal modes for which the ratio (B/a) 
does not change with time. 

The general solution can be expressed as a linear combination of 
these modes that satisfies appropriate initial conditions. Two factors 
complicate its specification. In the first place, the isotropic density 

. iu) t Im t 

die e^'"M(<Qo>cx+ +<Qj^> e+)e + ( < Qq > a_ + <Q^> 8_>e ' ] 

— 00 

(42) 

which is what experiments measure, should be described rather than the 
individual components f^ and In the second place, the initial aniso- 

tropy must be more carefully treated in focused transport, where it has 
an important effect upon the evolution of the density, than in purely 
diffusive transport, where it has a minimal effect. The solutions given 
below describe how the density depends upon space and time after the 
injection of a localized pulse with a finite initial velocity. In this 
situation, the injection velocity specifies the initial anisotropy in a 
physically meaningful way. 

The initial conditions that correspond to such an injection at t = 0, 
z = 0, of a pulse moving in the +z direction with a velocity are first 

Fq{z,0} = Nq 6{z) = J dx e , (43) 

in which 5{z) is the Dirac delta function, and second 
F^{z,t) “ e"^^^ F{z - V^t} 

in which the exponential factor is necessary for proper normalization and 
which is equivalent to 
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dF 

3t^“ “ lJ^O • ( 4^) 

The linear combination of modes that matches these conditions must satisfy 
two equations involving the amplitudes a^, a , and B , 

[<Qo>“++<Qi>e+] + KQo>«. +<Qi>e_] = Nq/ 2 it , (45) 


and 

+ V^(k - i/L)] [ <Qo>a^. +<Q^>6^I 

+ + V^(ic - i/L)] [<pQ>a_ +<Qj^>6_] = 0 , (46) 

which correspond respectively to the first and second conditions. When 
the solution of these equations is substituted in equation (42), is 
given by 


lu t 


lu.t 


\ - i/L)]e ‘ - [U)_ + V,(ic - i/L)]e 


0 2it 


d< e 


U)^ - U) 


(47) 


which is similar in form to equation (11-47). 

With the aid of equation (36), this expression can be rewritten as 


N. 


0 2V^ 


exp{ - (t/o ) - K (z + V t)} + V^(K- - :- ) + 

* cLo i2L 


(48) 


+ ■'2''c I e - * V fe] . 

in which the propagation function P is defined by 

, 1 /•+i" exp{V^t(K ^ + s^)"'^} - exp(-V.t(K,^+ 


/> 2 ^ IM 

(Kj^ + s ) 


" V'l^V “ (z + V^t)"'^}, - V_t < z < + V^t, 

“0 , z < - V_t ; z > + V_j^t , (49) 

where s = - ix - tC 2 is the Laplace transform variable Invoked by Abramowltz 
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and Stegun (1964, eq. [29,3.92]). Near the origin, P behaves as a 
modified Bessel function, but it jumps to zero at two discontinuities ,one 
of which moves in the +z direction with velocity 

V. - V. - V , (50) 

+ * C 

while the other moves in the -z direction with velocity 

7 - V . + V . (51) 

_ * c 

At the discontinuities, the numerical value of P remains constant at 

unity. Between them, it grows monotonically with time to form a bell 

shaped spatial profile whose peak moves in the -z direction with velocity 

V . Thus, the phenomena predicted by equation (48) are qualitatively 
c 

similar to those discussed in Paper Il-SIII. More specifically, the 
spatial and temporal derivatives of the step discontinuities in P give 
rise to S functions which represent two localized pulses moving in 
opposite directions. Initially, these coherent disturbances contain all 
of the particles injected, but as trajectories are scattered, the number 
of particles in an extended wake, which is spread continuously between 
the pulses, grows larger. This wake, which arises from the continuous 
portion of P, develops into a moving Gaussian analogous to the familiar 
diffusive profile. 

Because of the artifact mentioned above, equation (47) is not 
properly normalized to a constant total number of particles. Thus, when 
the Integration specified by equation (5) is performed, the Integral 

over z gives a delta function in tc, 5{i/L}, which Immediately leads to 

, lu)_{i/L>t . . , i(*M-{i/L}t 

0 ) {1/L}e “ (D li/L>e 

-± = . (52) 

(i) {i/L} - (I) {i/L} 

▼ ^ 


r 


FQ{z,t} dz = N 
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Because equation (36) specifies that the frequencies u>^{i/L} and b>_{l/L) 
are both positive Imaginary numbers, this expression represents a sum of 
two exponentials decaying with different time constants. Similar sums 
arise in the radioactive decay of genetically related nuclides and In the 
fragmentation of cosmic-rays, but In contrast to these situations where 
the initial number of secondary particles is zero, the coefficients in 
equation (52) are such that the initial slope is zero. Consequently, 
there is a brief period after injection during which the normalization 
Integral is virtually constant. Then, it decays at an exponential rate 
corresponding to ui^{i/L} which is the smaller frequency. Figure 5 shows 
how the ratio w_{l/L} of this rate to the rate which characterizes 

the overall evolution of the distribution function, depends upon (V/AL). 

In all cases, the normalization failure is insignificant in the weak 

4 

focusing limit where this ratio varies as (V/AL) . In the case of 
isotropic scattering, q = 1, the normalization failure becomes intolerable 
in the strong focusing limit where the ratio approaches unity. In the 
case of anisotropic scattering, on the other hand, the normalization 
failure is not very important, for the maximum relative decay rates are 
only 8% and 5% at q = 1.5 and q = 1.9, respectively. However, even these 
small deviations from proper normalization can lead to significant effects 
at long times after the injection. Thus when t >> l/ii)_{i/L} and also when the 
scattering is isotropic and focusing is intense, the solutions given here 
are not accurate. Under these circumstances, higher-order components, 
such as I 2 and f^, should be taken into account. 

The average density defined as an unweighted integral of 

Fq over z. In this integration of equation (47), the factor exp{z/L) does 
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not appear, 


<^ o > 


and the delta function becomes d{0} which leads to 

- f FQ{2,t)d2 * Nq - Njj(V^Oj^/L)(l - e'*^^°l) . 

M0O 


(53) 


The average density is proportional to the numbe ' of particles injected, 
but it also depends upon initial conditions and time. The latter depend- 
ences are not surprising, because the local density depends, as a result 
of guiding field convergence, not only upon the number of particles 
present but also upon where they are located. For isotropic Injection, 

= 0, the average density is the same as if focusing were absent. This 
average density also occurs just after injection with a finite velocity, 
but the transient component decays with time constant to give at 
equilibrium an enhanced density for V. < 0 and a reduced density for V, > 0. 
An interesting example of this behavior is the case = (L/Oj^) in which 
the average density decays to zero. 

In Paper II, the density f^ was expressed as the following weighted 
sum of two elementary disturbances f^^ and f^ associated, respectively, 
with coherent pulses moving in the +z and -z directions: 


fQ{z,t,vJ= h[i + (v^/VQ^)]fo ■*' ~ ^V'^Ol^^V 


where 22 2 h 


(z + V t)^ 


'01 
(V^jit - z) 




1 - 


• N(je 


-till 


■I 


6(z + Vjj^t} + 


1q ((l/2T^Vp^)(V^j^ t - z ) ) 




+ V„jt) 


4t V 
1 01 


1 


(54) 
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which is equation (11-51) repeated to facilitate comparison. The 
analogous result obtained from equations (48) and (49) is 

FQ{a,t,V^} » hli + (V^ + VJ/VJFq'*' + H(1 - (V^ + V^)/V^]Fq“ 


where 


- Nq exp{ - (t/o*) - <2^* 

3 1 1 <i (z + V_t)’^ y, 

^ ^ ^ ^ Vz ■*' Z (y j. _ V^^JJ ’ 

Fq~ - Nq exp{ - ^ U 6{z + V_t} 

* 2^; I ^ * V2 * ^<’'2 - L»'o'>'> * - ;T77^ V’'»U • 


and where y is the argument of the Bessel function that appears in 
equation (49) . The rectilinear transport described by the first of 
these expressions differs from the focused transport described by the 
second in the following respects; 

(1) The weighting factors H(1 - (V^ + V^)/V^] are such that the 
contribution of F„ vanishes when V, = V. while that of F„^ vanishes when 

U T + U 

*4“ “ 

= - V . The corresponding disappearance of f^ and occurs at equal 
positive and negative velocities, V, = + V*.. However, in all cases, if 

I Oi 

the injection -elocity coincides with the velocity of either coherent 
pulse, then Liie other pulse is absent. 

(2) In addition to a temporal dependence similar to that in 
equation (54), the exponential factor in equation (55) also embodies a 
spatial dependence such that the pulse moving in the -z direction is 
enhanced relative to the one moving in the +z direction. 

(3) In contrast to and f^” which are independent of the 
injection velocity, the elementary disturbances F^ and Fq depend 


(55) 

I 
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explicitly upon through the factor multiplying the 1^ component of 
the wake. 


(4) Xn the arguments of the Bessel functions and in the factor 
multiplying the component, the parameter (1/2 t^V^^) that appears in 
equation (54) is replaced by in equation (55). 

Because of its complexity, the implications of equation (55) are 
best discussed, as is done below, in terms of limiting cases for which 
a simplification occurs. To give perspective on these illustrations, 
figure 6 shows how the velocities of the coherent pulses depend upon 
(V/AL) . When the scattering is anisotropic, as it is at q 1.5 and 
q ■ 1.9, V is slightly larger than V , and in the strong focusing 

limit, both velocities Increase slowly with (V/AL). However, when q = 1, 
the coherent velocities diverge, for V increases as before while 
continues to decrease with (V/AL). Perturbation theory led to the 
identity “ (1/2L). The approximate validity of this relationship is 
demonstrated, for q = 1.5 and q = 1.9, in figure 7 where the product 
lies within a few percent of unity over a wide range of (V/AL) . 

The dotted curve, which refers to the co-ordinate scale at the right, 
shows that the relationship is also valid within + 20% for q * 1.0 
provided that (V/AL) < 5. 

The Gaussian limit of equation (54) applies when V^^^t >> [z| and 

t/2j^ » 1. Under these circumstances, in which the Bessel functions 

h 

approach their asymptotic form exp{y}/(2iry) the square root appearing 
in their arguments can be expanded with the aid of the binomial theorem 
to yield 


f^Ca.t) = 


Fp exp{- z /4Dt) 
2(TiDt)'^ 


(56) 
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which Is the tomlllar Green's function for one dlnenslonal diffusion from 
impulsive injection at t " 0 localised at s = 0. and which involves the 
coefficient of diffusion 0 = (See e,. (11-43] for a more direct 

derivation of this expression.) Similarly, when V^t >> Is + V^t| and 
u t » 1, and if the second terms in the asymptotic expansions of 1^ 
and l^ (Abramowits and Stegun, 1964, e,. (9.7,11) are also taken Into 
account « equation (55) reduces to 


^ "f" 

F ■ N_C exp{— + 

*^0 0 Til t ^ 



(z + V^^t ~ X) 


exp 

4Lj,t 






2(iTD^t) 

-■ 


(57) 


where 




2k. 


u. 


(58) 


'1 '"01 

18 an Important new parameter that I call the coefficient of focused 
diffusion, where 

19 the velocity with which the peak of the bracketed Gaussian moves in 


(59) 


the -z direction, and where 

1 ' * 


(V* - 


(60) 


T„ Ao.V^U 

tf 1 * 01 

describes an exponential growth that does not occur in the absence of 
focusing. Similarly, the factor 

2 


r 


1 + 






( 61 ) 


depends upon and deviates from the value of unity that it has in 
equation (56) . 
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The length X, given by 
a^V,(V^ + VJ 


c u 


(62) 


01 


represents a small displacement of the z co-ordinate which arises because 
Injection occurs, in effect, over the distance covered by the 

initial pulse before it decays rather than strictly at z =0. The time 


T^, given by 


!i!^ 

" A U 


01 




'V + V, 
c ■ 

cv^ 


.2 1 


(63) 


characterizes a small correction which depends upon the injection velocity 
but which becomes negligible as t increases. These two corrections arise 
from the second terms mentioned above. Although they do have a minor 
effect upon rectilinear transport, they do not appear in the standard 
Gaussian approximation where isotropic injection is implicitly assumed. 
They were taken into account here in an attempt to improve the accuracy of 
the Gaussian representation, but for many purposes they can be neglected. 

Equations (57) -(60) have three Implications which mean that focused 
diffusion is strikingly different from ordinary diffusion. In the first 
place, the point of maximum density, which occurs at z = X - V^t where 
the argument of the Gaussian function is zero, moves in the -z direction 
with velocity In ordinary diffusion, the Gaussian remains centered 

at z = X. In the second place, the coefficient is substantially 
larger than D whenever (V/AL) is large. In the third place, th-=i factor 
C exp{t/t^} multiplying the Gaussian not only depends upon the injection 
velocity V^, but also increases exponentially with time. 
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provide a reference to ordinary diffusion against wbi.ch these 
implications can be compared, figure 8a, at the left, shows profiles of 
density vs z, for q - 1.5, at the instant t = (10/A) « following an 
injection with =■ + The wake calculated from equation (54) , which 

is rigorously exact and which is shown as a solid line, drops dlscontinuously 
to zero at z » ± ^01^* Because the velocities at injection are collimated 
in the +2 direction, the wake peaks at a point, z * 

right of the point of injection. Following a suggestion made in 
Paper II, the coherent disturbance (dotted curve) is represented here by 
a Gaussian whose width (D^t)^ is characterized by a coefficient of 
dispersion D* = (D/20). Within the region where the density in the wake 
is finite, the dashed line, which represents the Gaussian predicted by 
equation (57) for L = », provides a fairly accurate description of the 
actual wake. Near the maximum, the two curves are nearly coincident, but 
the Gaussian is about 25% too low at z = + V^^t where the asymptotic 
expressions for the Bessel functions are not accurate. For iz| > V^j^t, 
the Gaussian gives a finite density which misrepresents the actual value 
of zero. Nevertheless, the contribution of this incorrect prediction to 
the total area under the dashed curve compensates for the underestimated 
density in the range UI < V^j^t in such a way that the area under the 
Gaussian is the same as that under the exact profile. Thus, in the case 
of ordinary diffusion, the Gaussian approximation is correctly normalized. 
For injection toward the left with the profile is the mirror 

image obtained by reflecting around z = 0 the profile for 

In figure 8b, profiles are shown for focused diffusion with (V/AL) « 
1.0, q - 1.5, t - (10/A) and = + V_^ = 0.549V. Because these profiles. 
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at the center, refer to the same time and to nearly the same positive 
injection velocity as those for ordinary diffusion, at the left where 
* 0.563V, the qualitative differences that are Immediately 
apparent must be attributed to focusing. Even though the injection is 
toward the +z direction, the wake is largest near the discontinuity at 
z ■ - V t. Near the other discontinuity at z = + V t, the density in 
the wake is relatively small. Consequently, the coherent pulse, which 
is represented as before by a Gaussian plotted as a dotted line, seems 
insignificant even though it contains 53% of the injected particles. 
Because the velocity = 1.081V is greater than V_ = 0.631V, the 

Gaussian peak lies far to the left of V t. Consequently, the density 
predicted by equation (57) (dashed line, for which = 1.07(V^/A) - 
1.34D) decreases monotonically with increasing z from its maximum 
value at - V_t. This decrease does not reproduce the exact prediction 
of equation (55) (solid line) which exhibits a maximum just to the 
right of - V_t. Nevertheless, the largest difference between solid and 
dashed curves is only 6%. It Is apparent in figure 8b, where the 
vertical scale is expanded by a factor of 3 over that in figure 8a, 
that the density in the wake for focused diffusion with positive inject- 
ion velocities is much smaller than that for ordinary diffusion. 

In figure 8c, which refers to the same conditions as figure 8b 
except that the injection velocity is negative, V^ = - V , the exact 
profile (solid line) decreases monotonically from its maximum at 
z ~ - V_t. The wake, which is plotted on a vertical scale reduced 
relative to that of figure 8b by a factor of 10 and relative to that 
of figure 8a by a factor of 3.33, is much larger than the wake associated 
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with ordinary diffusion or with focused diffusion after injection 
toward diverging guiding fields. The coherent pulse (dotted line), 
which is 313 times larger than the one in figure 8b, is a prominent 
feature of the density profile, but it contains only .12% of the injected 
particles. The Gaussian approximation (dashed line) peaks well to the 
left of z ■= - V t . It approximates the solid profile with a 

maximum deviation of 22% . 

Because the bracketed expression in equation (57) has the form of 
a normalized Gaussian, the expression multiplying it represents the 
total area under the approximate profile from z=-‘»toz=+“. 
Obviously, from figures 8b and 8c, this area has no straightforward 
relationship to the area under the exact profile, for most of it lies 
near the peak in a region where the density actually is zero. Con- 
sequently, the expression multiplying the brackets should be regarded 
as one which gives an accurate approximation and not as a normalization 

parameter. 

The tendency of the Gaussian peak to outrun the discontinuity 
moving in the -z direction must disappear when focusing is very weak, 
for the unfocused Gaussian is stationary. This disappearance is 
illustrated in figure 6 where the dotted 11;.. ...presenting crosses 
the solid line representing V_ at a value of (V/AL) below which focused 
diffusion is governed by the extreme weak focusing limit discussed in 
the next paragraph. Evidently, this value becomes small as q - 2. 

To Interpret solar particle events, many authors have assumed that 
Interplanetary propagation is governed by the diffusion equation for 



31 


spherical geometry. In the present context, this equation corresponds 
to the result of substituting Into equation (6) a first-order approxi- 
mation to the flux, S = - D(3Fg/3z), (see Paper III) to yield a diffusion 
equation 



3t L dz 



( 64 ) 


With the aid of the transformation 
w » 2 + (D/L)t , 

equation (64) can be put In a form. 


!!c 

at 


2 


3w 


( 65 ) 


whose solution, 

p = U exp{ - w^/4Dt> _ exp{ - [z -f (D/L) t] ^/4Dt) 

° ° 2(irDt)^ ° 2(iTDt)''^ 


( 66 ) 


is a Gaussian similar to the one in equation (57) except that D appears 
in the place of Equation (66) does not take into account the 

quantitative difference between these coefficients, it does not Include 
the exponential growth that occurs when focusing is intense, and it does 
not describe either the prominent coherent pulses or the pronounced 
dependence upon •f- . tion velocity that characterize focused diffusion. 

Because these considerations also apply to the diffusion equation for 
spherical geometry, existing work on interplanetary diffusion should be 
re-examined from a point of view which correctly incorporates the effect 
of focusing. 

The limit <j^V^t 0 is a simple case that illustrates the transition 
from focused diffusion to supercoherent propagation but that does not 
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depend upon the detelled properties of Bessel functions. In this 
special case, where Iq = 1 and “ 0, e<{uation (55) reduces to 

= Nj e«p{ - (t/o,) - <^(z + V^t)) jtiU + (V^ + V^)/V,)6{z - V^t) 

* Ml - (V^ + V^)/V,1«P + V_t) + 2^ [^ - 

in which the identity ® (1/2L) has been invoked. Here, the coherent 
6 functions are weighted as before, and the wake spread between them 
depends upon z and t only through the exponential multiplying the 
brackets. In the diffusive regime, the Gaussian form reached by the 
wake is virtually independent of conditions at injection, but in the 
relatively strongly focused regime exemplified by equation (67) , the 
wake depends critically upon V.^. In particular, if and 

L ^ a. (V. - 2V ) , the wake is completely absent , and there is only 
one coherent pulse. It moves toward regions of reduced guiding field 
while its an^litude decreases as exp{ - (z/L)}. Because the exponential 
factor appearing in equation (5) compensates for this decrease in 
density, the total number of particles in the pulse is constant. On 
the other hand, if = - V , there is one coherent pulse moving with 
constant amplitude into stronger guiding fields. There is also a 
substantial wake whose density decreases exponentially with distance 
away from the pulse. In this situation, according to equation (5) , the 
number of particles in the pulse decreases exponentially with time, but 
this decrease is accompanied by a growth of the number in the wake 
such that the total remains constant. 

As was discussed in Paper II-5III, the first of these examples 
embodies the most pronounced alignment in the +z direction of particle 
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velocities at injection that can be adequately treated in a description 
which invokes only the two components and In this situation » 

the wake is Insignificant and nearly all of the injected particles 
remain in a coherent pulse whose amplitude decreases rapidly as a result 
of the geometrical divergence of field lines but very slowly as a result 
of scattering. Qualitatively, this stability of the number of particles 
within a bunch whose velocities are collimated along the guiding field 
occurs because the rate at which they are realigned by adiabatic 
focusing exceeds the rate at which they are scattered. In the second 
example, on the other hand, the same effect works in the opposite 
direction, for focusing aids scattering by rapidly removing particles 
from a bunch whose velocities are aligned in the -z direction. In spite 
of this rapid reduction in the number of particles, the density within 
the bunch remains stable, because particles moving in this direction 
converge together laterally along with the guiding lines of force. 
Because particles removed from the pulses constitute the source of the 
wake, the above considerations also explain why the wake associated 
with an injection toward stronger guiding fields is more pronounced 
than the one associated with an injection toward weaker fields. 
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V. THE SUPERCOHERENT MODE 

The discussion In §IV outlined the changes which occur as the intensity 
of focusing is increased while other conditions are unchanged. In the weak 
focusing limit described by equation (66) , where the scattering length (V/A) 
is much smaller than L, the Gaussian profile of ordinary diffusion drifts 
slowly into stronger guiding fields. In focused diffusion described by 
equation (57), where (V/A) « L, the drift velocity exceeds the particle 
velocity, the diffusive wake, which is no longer a bell shaped Gaussian, is 
dependent upon conditions at injection, and the coherent disturbances are 
very prominent. In the strongly focused regime described by equation (67), 
where (V/A) j L, coherent effects dominate and the wake is Insignificant. 

In the supercoherent mode to be discussed in this section, which occurs 
when (V/A) » L, the tendency of focusing to enhance coherent transport and 
to suppress diffusive transport reaches a limit in which the wake is 
completely absent and particle propagation is coherent. The word super- 
coherent is appropriate here because transport phenomena in this regime are 
analogous to those in the superfluid and superconductive states. Unlike 
these states, the supercoherent mode does not appear at a discontinuous phase 
transition. However, the supercoherent transition does occur very abruptly 
with the disappearance of the velocities in figures 3b and 3c which embody 
the coupling between coefficients of opposite parity that leads to diffusive 
effects. 

Because the three negative velocities in figures 4b and 4c display almost 
the same pattern as the three positive velocities, it might be expected that 
a second supercoherent mode would propagate in the direction opposite to the 
one that actually does propagate in the +z direction. In fact, the second 
mode decays rapidly and leaves its particles in an extended wake. The 
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reason for this decay can be seen in equation (25) where the expression 
within square brackets that multiplies (3£^/9z) must vanish according to 
equation (26) . But the terms proportional to and are small when 
(V/AL) Is large, while those proportional to and 0^2 are large and 
negative. Consequently, because the truncated form of this expression takes 
on negative values instead of the required value » 0, the flux equation 
la not even approximately satisfied when f^ and are finite. Thus, a 

mode that propagates coherently in the -z direction can be constructed as 
a linear combination of and Q 2 , but it can not be properly normalized. 

In this situation, which involves strongly focused disturbances that move 
In the -z direction, a perturbation approach analogous to that employed in 
Paper III is more appropriate than the method of eigenfunctions. Such an 
approach cannot be pursued here, but its qualitative effect can be judged 
In Figure 8c where dispersive effects would mix together the triangular 
wake and the coherent pulse to give a broad disturbance propagating quasi- 
coherently in the -z direction with an effective velocity of I call 

this the pseudodif fusive mode. 

Only a few of the positive terms that appear inside the bracketed 
expressions multiplying (af^^/Sz) and Of2/3z) in equation (25) are needed 
to approximate and V^. Consequently, an accurately normalized super- 
coherent mode can be constructed as a linear combination of Qj^ and Q^. 

Because these eigenfunctions are not coupled to and Q 2 , only the two 
odd transport equations are required to describe this mode. This simplifi- 
cation is similar to the one that appeared in Paper II-5IV where the purely 
coherent modes that occur when q > 2 were discussed in' terms of two transport 
equations. However, these modes, which arise because pitch angle scattering is 
very weak near y = 0, are physic.ally different from the supercoherent mode. 



36 


which occurs because the tendency of focusing to align particle velocities 

overcomes the tendency of scattering to make them isotropic. 

The equations which describe the supercoherent mode were obtained from 

equations (16) and (18) by neglecting the small gradient terms in 

^23 retaining the large ones in and They are, 
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Because their form is very similar to that of equations (29) and (30), 
there is no need to derive in detail the solutions discussed below, for 
each step corresponds exactly to a step taken in §IV to derive solutions 
for focused diffusion. Thus, the frequencies and can be obtained by 
substituting into equation (36) the following expressions! 

(70) 


°3 
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- * “n' • 
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Because the parameters that specify and oj are the same as those that 
appear in equation (55), these redefinitions determine completely the 
exact solutions that apply in the supercoherent regime. 

If these frequencies are substituted in equation (52) , it predicts truly 
minuscule deviations from proper normalization. In contrast, the frequencies 
that apply to the mode that propagates in the -z direction, which follow 
from two transport equations obtained from equations (15) and (17), 
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lead to an extremely rapid decay of the normalization Integral. This 
analysis establishes the posit ive**veloclty supercoherent mode as a unique 
and fundamental feature of strongly focused transport. 

In spite of the mathematical identity embodied in equations (70) -(74), 
supercoherent propagation is different from focused diffusion. In the 
latter regime, where V is small and V. is large, the bell shaped profile 
of the propagation function drifts slowly in the -z direction while the 
discontinuities move rapidly in opposite directions. This configuration 
leads to the diffusive evolution of a wake spread between two coherent 
pulses. In the former region, where is large and is small, both 
discontinuities of the propagation function move in the +z direction, and 
the bell shaped profile spread between them also moves in the same direction 
with an intermediate velocity. This configuration leads to the unidirectional 
propagation of localized disturbances. The velocity V of the leading 
discontinuity ahead of which there are no particles is given by 


^ “l3> * 


( 77 ) 


while the velocity V, of the trailing discontinuity behind which there are 

b 


no particles is given by 


*=<"11 ^ “ u ' - 


( 78 ) 


If the injection velocity is positive and if it lies between and V^, 

the density profile is given by equation (55). This rigorous solution 


includes a coherent disturbance associated with each discontinuity, but 
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these higher-order pulses, which are analogous to those discussed in 

Paper II-IIV, rapidly become insignificant compared to the continuous 

profile. The breadth of this profile can be determined from figure 9 

where the velocities of its boundaries, V and V. , are plotted against 

a D 

(V/AL) for q = 1.5 and q - 1.9. If > V^, ev-,uatlon (55) does not hold, 
because eigenfunctions above must be included in a proper description 
of the strong anisotropies implied by this initial condition. However, 
these anisotropies decay rapidly to leave a situation similar to the one 
that follows the injection of a pulse with velocity V . Because the same 

a 

consideration also applies to the strongly anisotropic Injection of solar 
particles, profiles calculated with are most appropriate for 

comparison with observations. If < V^, the supercoherent disturbance, 
which is very similar to the one that follows injection with Vj, = is 
accompanied by a broad pseudodif fusive wake moving in the -z direction. 

The proper description of this disturbance also involves eigenfunctions 
above Q^, but negative or small positive injection velocities do not 
correspond to solar injection. Thus, for virtually all plausible initial 
conditions, the supercoherent profile can be accurately described by 
substituting the parameters defined by equations (70)-(74) in equation (55) 
The Gaussian limit of this solution, which applies when V^t » jz - V^t 
is 


Pq = Nq C exp{ 



where 
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is the coefficient of supercoherent dispersion which plays much the same 
role as the coefficient of dispersion defined in Paper II-5IV, where 


- V - — V* “ H(1 

c * 


B„-> “U * - U^ 3 > "33 


( 81 ) 


is the supercoherent velocity with which the peak of the Gaussian moves in 


the -fz direction, and where 

1_ i_ - "l3>' 

'S “1 ”3 ”1 “X 3 


(82) 


describes an exponential decay of the amplitude. The factor 




1 + VI ~ ''3/ 




J . MBj/ - U33') - >,V, I - U33 . 


V.Bi3 


(83) 


embodies the dependence of the supercoherent Gaussian upon injection velocity. 
The offset is given by 


X 



V*(V^ + V^) 

U13C 


(84) 


The correction characterized by Xj. is negligible. 

In figure 9a, the dotted line which represents for q = 1.5 almost 
coincides with the dashed line which represents V^. In figure 9b, where 
q = 1.9, this coincidence is r.''.arly exact, but, for clarity, only the dotted 
line is shown. This means that the supercoherent Gaussian moves with virtually 
the same velocity as the coherent delta fur.ctlon of focused diffusion. In 
Paper 1I-5IV, the coherent Gaussians, which had the same velocity as the 
coherent delta functions, embodied an improved representation of coherent 
disturbances in which dispersion was included. Similarly, it appears here 
that the consideration of higher-order eigenfunctions leads to an improved 
representation which embodies the dispersive evolution of the supercoherent 
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dlscurbance. To carry this line of reasoning one step further, if the super- 
coherent pulse is equivalent to the delta function in equation (67), then the 
wake predicted there should give a reasonable estimate of the small diffusive 
wake that remains in the supercoherent regime. 

These points are llixostrated in figure 10 where a supercoherent pulse 
for q = 1.5, (V/AL) “ 5, => » .8V and t * (2/A) is shown by a solid line, 

which represents the exact solution, and by a dotted line, which represents 
the Gaussian approximation. Associated with this pulse is a wake predicted 
by equation (67) which is also shown as a solid line, because <2 is negative, 
exp{ - <2 increases with z. Because of this weighting, the supercoherent 
velocity is only slightly less than the velocity of the leading dis- 
continuity, and the profile drops to zero just in front of its peak. Con- 

l- 

sequently, the Gaussian width (Dgt) ^ overestimates the actual width. However, 
this complication does not occur when t The same weighting effect^which 

puts the Gaussian peak in a region where the asymptotic representation of the 
Bessel functions is not accurate, underlies the 'v 15% difference by which the 
exact and approximate peaks are separated in figure 10. 

In figure 10, the dashed wake and Gaussian pulse are those of figure 8b with 
their horizontal scale transformed In such a way that this profile for focused 
diffusion corresponds to the same time and to the same value of L as the 
supercoherent profile, but to scattering 5 times more intense. This change in 
A, which goes across the supercoherent transition, leads to a dramatic increase 
in the magnitude of the wake relative to that of the pulse. However, the 
dotted supercoherent Gaussian has about the same shape as the dashed coherent 
one. More quantitatively, the effect of Increasing (V/AL) from 1 to 5 is to 
reduce the coefficient of dispersion from a virtually unfocused value 
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- .0405 (V^/A) to the supercoherent value Dg = .0123 {V /A). This 
reduction by a factor of 3.3 is slightly overcompensated by the factor of 5 
decrease in A, but the supercoherent Gaussian is not perceptably wider because 
the width, at a given time, has a v^ak square-root dependence upon the 
coefficient of dispersion. 

The Gaussian supercoherent pulse embodies an equilibrium pitch angle 
distribution in which the opposing effects of scattering and focusing balance. 
Because of the collimation produced by focusing, stochastic variations in the 
pitch angle of an individual particle average to give a net velocity which 
is finite and approximately equal to that of the other particles in the bunch. 
Statistical fluctuations in this averaging of random velocities give rise to 
dispersion. This section has put this physical picture on a rigorous basis. 
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VI. INTERPLANETARY PROPAGATION OF ENERGETIC PARTICLES 

Tweaty years ago, Meyer, Parker and Simpson (1956) concluded chat 
diffusion could explain the temporal profile of solar cosmic-ray intensity 
on 23 February 1956. Subsequent investigations have confirmed the basically 
diffusive nature of particle propagation in interplanetary space, but they 
have also uncovered many effects that can not be understood in terms of 
pure diffusion. The objective of this section is to show that several of 
these unexplained features arise as natural consequences of ad-abatlc 
focusing in the spiral Interplanetary field. Because the present theory 
takes Into account neither perpendicular diffusion nor convenctlon, it 

would be premature to attempt a quantitative comparison of observed and 
predicted solar event profiles. Instead, the discussion that follows gives 
a qualitative interpretation of prompt events in which these relatively slow 
processes play a minor role compared to rapid propagation along field lines 
that trace out a reasonably direct connection between the Earth and a flare 
on the western limb of the Sun. Focusing must have important effects in 
cosmic-ray modulation, but this steady state phenomenon, in which convection 
plays a crucial role, also lies beyond the scope of the present discussion. 

One of the most striking and least understood aspects of interplanetary 
physics is the "scatter- free** propagation of kilovolt flare electrons in 
which an impulsive burst of particles, usu;>lly followed by a slowly decaying 
tall, arrives at Earth with an average velocity parallel to the field of 
^0.8V (Lin 1974). From observations of Type III radio noise generated by 
this bunch of electrons, Lin, Evans, and Fainberg (1973) have demonstrated 
that it travels »1.2 AU along a spiral field line. Although this phenomenon 
corresponds to the coherent mode discussed in Paper II, there are two dis- 
crepancies which make untenable an interpretation based upon this agreement. 
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In the first place, the predicted velocity is closer to 50% of the particle 
velocity than it is to 80%. In the second place, developments in scattering 
theory subsequent to Paper II indicate that scattering near u = 0 is not 
weak enough to allow the purely coherent mode to persist (Jones, Kaiser and 
Birmingham 1973, Volk 1973, Owens 1974). These discrepancies do not apply 
to the supetcoherent mode, for it occurs even in the presence of scattering 
at w - 0,’ and its velocity is close to the observed one. Thus, we can 
interpret "scatter-free" events as supercoherent bunches which propagate with 
very little dispersion in the strongly diverging fields near the sun. In 
this situation, the parameter (V/AL) is not constant, for the scale length 
of the interplanetary field is approximately (r/2) where r is the distance 
to the sun. Moreover, the radial dependence of (V/A) undoubtedly leads to 
additional variations of the focusing parameter. If (V/AL) lies above the 
supercoherent transition, these variations merely cause the supercoherent 
velocity to change without affecting the basic nature of the mode. This 
decrease in parallel velocity corresponds exactly to the "deceleration of 
the Type III exciter reported by Evans, Fainberg and Stone (1973J. However, 
it occurs because the pitch angle distribution of particles with a given speed 
becomes broader as focusing decreases, and not because the speed of an 
individual particle changes. Thus, the observed transit time, which represents 
an average over a gradual decrease in velocity, corresponds to a velocity 
larger than the local one. However, this enhancement should be small, 
because regions where the velocity is slow are more heavily weighted than 
are those where it is fast. Most supercoherent events involve electrons, 
but they can also involve protons. An example is the event on 24 March 1966 
(McCracken, Rao and Bukata 1967). 

Within the same framework, the properties of two other types of solar 
burst, reviewed by Kundu (1965), can be understood. The U-type bursts, in 
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which an ascending frequency branch appears after the descending branch of 
a Type III burst, are interpreted in terms of closed magnetic lines along 
which the electron bunch returns toward the sun after moving up from the 
flare site. The "shortness" of the ascending branch compared to the 
descending one can be explained In terms of the greater dispersion of a bunch 
moving pseudodlf fusively down into converging lines relative to that of a 
bunch moving supercoherently upward. The V-type bursts, in which a Type III 
burst Is followed by a brief period of broad-band continuum radiation, can 
be explained as synchrotron radiation from a nearly isotropic cloud of 
electrons which forms after the bunch passes through the supercoherent 
transition. This phenomenon is very similar to those discussed below at 
greater length. 

The chief objection to this interpretation is that the values of (V/AL) 

calculated for electrons from observations of interplanetary field fluctuations 

lie below the supercoherent transition. Thus, in the example discussed in 

Paper II and presented there as figure 6, where 38 KeV electrons were 

-4 2 -1 

scattered by fluctuations for which P = 6.3 x 10 y H 2 at f„ = 0.5 Hz 

XX 0 

and for which q = 1.9, the focusing parameter is (2V/Ar) = 0.022 which is 
well below the value (V/AL) 0,5 required, in figure 3c, for supercoherent 
propagation. Within the framework of current theories which predict that 
pitch angle scattering increases with decreasing rigidity, this objcctfon .an 
not be overcome. However, there are compelling intuitive reasons to expect 
a regime where low rigidity particles follow adiabatically the stochastic wander 
of magnetic lines of force with very little scattering relative to the local 
field direction. The correspondence between the "scatter-free" and super- 
coherent effects can be regarded as empirical proof of the existence of such 
a regime where scattering decreases with decreasing rigidity as it merges 
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into the Alfvenlc idealization. Consequently, more theoretical effort should 
be devoted to exploring the low rigidity limit of pitch angle scattering. 

Less effort should be devoted to analyzing the detailed dependence of ^ upon 
p, for most of the special significance of weak scattering near p = 0 
disappears in the present context where focusing moves particles rapidly 
through this region. 

There. is a localized region in which the nature of interplanetary 

propagation changes from supercoherent to diffusive, for the focusing parameter 

must eventually pass through the supercoherent transition by virtue of its 

inverse dependence upon r. In this situation, shown schematically in 

figure 11, solar particles propagate supercoherently to a fairly abrupt 

transition (wiggly line) beyond which focused diffusion occurs. In effect, 

the injection of particles into this region of focused diffusion is highly 

anisotropic and occurs far from the sun. Thus, because V, = ss V , the 

t § + 

+ • 

density profile can be described by the function given by equation (55) 
with the origin, z = 0 and t = 0, chosen to be the place and time at which 
the supercoherent pulse hits the transition. 

A well known prediction of the telegrapher's equation is that, within 2 
mean free paths of an impulsive injection, the temporal profile has an 
abrupt onset followed by a monotonic decay. At larger distances, the profile 
exhibits a relatively gradual Increase from onset to a broad maximum which 
is followed by a monotonic decrease. Similarly, there is a certain distance 
from the point of injection within which the maximum density predicted by 
Fq ^ 2 , t) occurs at onset and beyond which it occurs after onset. This 
distance is given by the positive root of the following quadratic equation: 



( 85 ) 
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Because Its coefficients are complicated, the implications of this equation 
are best illustrated by quoting the numerical result 2 ^ 35.8 (V/A) obtained 

for q - 1.5 and (V/AL) = 1, which is the case illustrated in figure 8. This 
distance is 6.3 times larger chan 2 mean free paths. Thus, a basic effect 
of focusing is to enlarge the spatial region within which temporal profiles 
have abrupt onsets. 

The region of focused diffusion can be divided by the dashed line in 
figure 11, which is located at a distance beyond the supercoherent 
transition, into two zones in which flare profiles are qualitatively 
different. The location of Earth relative to the wiggly and dashed dividing 
lines in figure 11 depends upon the intensity of interplanetary magnetic 
fluctuations and upon the velocity and rigidity of the particles being 
observed. Consequently, the nature of observed flare profiles is expected 
to show considerable variability from day to day, and, at a given time, it 
may not be the same for all particle species. 

In the zone between the dashed and wiggly lines, dispersive effects 
broaden the coherent pulse and smear out the discontinuous onset of the wake. 
Consequently, the temporal profile exhibits a fairly abrupt onset whose 
duration corresponds to the width of the coherent Gaussian. Then, the intensity 
decays monotonically from its maximum. Many authors have assumed that the 
interplanetary diffusion coefficient is the one which appears in the diffusive 
Gaussian that best fits the profile observed during a given event. But, if 
this procedure is applied to an abrupt-onset event, the coefficient of dis- 
persion is obtained Instead of the coefficient of diffusion. (See Paper II, 
eq. [94].) This misidentlfication explains why the coefficient of diffusion 
seems to be paradoxically small when the overall behavior of the profile 

corresponds to large mean free paths. 
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Abrupt-onset flare profiles are very commonly observed not only for the 
low-energy solar protons and electrons recorded by satellite monitors but also 
for the high energy nuclei sensed by ground level neutron monitors. Their 
hi^ probability of occurence is inconsistent with a purely diffusive picture, 
because the 2 mean free path zone where they are expected is narrow and lies 
close to the sun. Because focusing widens this zone and moves the point of 
injection outward toward Earth, the frequent occurence of abrupt-onset profiles 
is a natural feature of the configuration shown in figure 11. To describe in 
detail the strong initial anisotropies of such events (McCracken 1962, 

McCracken, Rao and Bukata 1967), it is necessary to evaluate the individual 
components f^ and fj^ whose separate behavior differs from that of the linear 
combination specified by equation (42) . This calculation can not be under- 
taken here, but the enhancement of coherent effects caused by focusing means 
that anisotropies persist for much longer times than they would in the absence 
of focusing. On many occasions, the initial anisotropy disappears suddenly 
about an hour after onset. This phenomenon can be understood from figure 11 
as the arrival of a disturbance propagating pseudodif fusively back toward 
the sun from distant regions of weak focusing where scattering finally succeeds 
In making the distribution function isotropic. 

Beyond the dashed line where maximum intensity occurs after onset, flare 
profiles resemble those of classical diffusion. If a diffusive Gaussian 
is blindly fitted to one of these profiles, which are described by equation (57), 
the best fit occurs, not for the coefficient of ordinary diffusion D, but 
Instead for the coefficient of focused diffusion Because > D, this 

mlsidentificatlon can lead to an overestimate of D. Moreover, in 

determining the distance to be invoked in the fitting procedure, both the 
distance to the supercoherent transition and the offset A should be subtracted 
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from the Sun to Earth distance. If these corrections are not subtracted, 
then D| is underestimated. 

In the late stages of many flare events, the decay takes on an exponential 
character which differs from the temporal power law predicted by diffusion 
but which can be explained by postulating a free-escape boundary far from the 
Sun (Meyer, Parker and Simpson 1956). Such an exponential decay appears as 
an inherent characteristic of focused diffusion. Its relaxation time t^, 
derived from equation (57), is given by 



This time decreases from its infinite unfocused value to approach 
in the strong focusing limit. Because typically observed relaxation times 
of many hours are much greater than any plausible value of c^, the rate of 
exponential decay must be controlled by weakly focused diffusion in the 
outer solar system. An accurate description of this regime would take into 
account the radial variation of the focusing parameter. This can not be done 
here, but it seems reasonable to expect that a region of uniform density would 
be set up in the inner solar system through the rapid equalization of density 
Inhomogeneities by coherent effects. Farther out, gradients would develop in 
a region of weakly focused diffusion through which particles escape at a 
rate leading to a slow exponential decay. 

The Jovian electron bursts recorded on Pioneer 10 result from the inward 
propagation of low rigidity particles along interplanetary fields (Chenette, 
et al. 1974, Teegarden, et al. 1974, Smith, et al. 1975). The supercoherent 
mode does not apply here, but quasi-coherent propagation into converging fields 
is predicted by equation (57), for relatively weak focusing, and by the 
pseudodif fusive idealization, for strong focusing. The brief duration and 
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10 hr modulation of the Jovian bursts can be explained in terms of these 
modes, but they are difficult to reconcile with pure diffusion. Outside 
Jupiter, electrons might propagate supercoherently . Here, the intensity of 
Jovian bursts would decay more rapidly with distance than it does inside, but 
their temporal fine structure would be better preserved. Observations of 
these effects could confirm the asymmetrical nature of focused transport. 

Several poorly understood interplanetary phenomena correspond to 
predicted features of focused transport. However, the current status of pitch 
angle scattering theory at low rigidities does not allow these features to be 
quantitatively related to the observed intensity of interplanetary magnetic 
fluctuations. In particular, a weak scattering regime at low rigidities is 
indicated not only by supercoherent electron events but also by the findings 
of Bryant, et al. (1965) in which profiles for proton events were dependent 
upon the distance travelled Vt but were independent of rigidity. Such 
behavior can occur only if A is a linear function of V alone, which means 
that (V/AL) is independent of both velocity and rigidity. Thus, the 
observations suggest that a broad region of rigidity-independent scattering 
lies between the low-rigidity regime, where scattering at a given velocity 
increases with rigidity, and the high-rigidity regime, where it decreases. 

The transport theory’ developed in this paper sidesteps these ambiguities and 
goes directly to the macroscopic description of density profiles in space 
and time. These profiles depend ultimately upon only three parameters, q 
which characterizes the anisotropy of scattering, A which characterizes 
the intensity of scattering, and L which characterizes the intensity of 
focusing. Until an improved description of pitch angle scattering comes forth, 
attempts to relate observed profiles to these parameters may be worthwhile. 
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VII, THE STRUCTURE OF EXTRAGALACTIC RADIO SOURCES 

Double radio sources, between which there Is usually found an 
optical galaxy, are among the largest and most energetic phenomena of 
astrophysics. In spite of the detailed knowledge of their structure 
" (Mackay 1971, Fomalont 1969) made available by recent advances in inter- 
ferometry (Ryle 1975), there is still no generally accepted explanation 
of these remarkable objects. The Interpretation put forth in this section, 
which adopts a widely held view that the radiating electrons gain their 
energy within the central galaxy, describes the symmetrical transport of 
these electrons to great distances from this source and the subsequent 
evolution of the clouds they form there, but it does not attempt to describe 
their acceleration. This Interpretation rests on the assumption that a large 
scale magnetic field, which threads through the galaxy, extends far into 
intergalactic space to form a diverging guiding field along which focused 
transport occurs. Here, the basic morphology of the double sources arises, 
in much the same way as in interplanetary propagation, when two bunches of 
electrons move rapidly out in opposite directions to supercoherent transitions 
where they form relatively long lasting clouds which constitute the actual 
radio sources. On the time scale implied by the large separation of the 
clouds from the central source, the Compton-synchrotron mechanism of energy 
loss plays an important role which will be mentioned below but which cannot 
be treated in detail. 

The chief objection to this interpretation is that the strong anisotropy 
of the supercoherent mode might be rapidly attenuated by the collective 
effects reviewed by Wentzel (1974). However, the conventional view, that 
these wave effects limit streaming velocities to the Alfven velocity , applies 
to steady state conditions. In the absence of reliable knolwedge of the inter- 
galactic medium and in view of the slow growth of waves, there is no reason to 
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believe chat scattering could be significantly enhanced during the brief 
period in which a supercoherent pulse sweeps over a given volume of the 
medium. In any case, as was documented above, the supercoherent mode can 
persist in the presence of scattering, including that generated by the 
particles themselves through collective effects. The interplanetary analogy 
may be relevant here, for Jovian bursts do persist while generating waves 
(Smith et al. 1975). 

Underlying the striking symmetry of double sources is the symmetry of 
the guiding field, which arises because Maxwell's equations guarantee that 
as many lines of force diverge out from a local condensation as converge 
into it. For example, if intergalactic currents are absent, the field would 
take on a dipolar character, and, if the central galaxy is surrounded by an 
expanding medium, the field would develop a radial pattern analogous to that 
of the interplanetary field. Within these bilateral configurations, 
electrons propagate and radiate in similar magnetic environments on 
opposite sides of the source. No matter how the electrons are accelerated, 
scattering within the central galaxy, which makes them isotropic there, 
ensures that the two bunches contain equal numbers of electrons. When these 
bunches hit the supercoherent transitions, particles are rapidly scattered 
into two clouds. Because these clouds contain equal numbers of isotropically 
distributed electrons and because their diffusive evolution is slight during 
the time required for light to travel between them, the two lobes of the 
radio source appear to have nearly the same luminosity regardless of the 
angle between their axis and the line of sight. The axis is perpendicular 
to the E vector of the polarized emission from the clouds. These predicted 
symmetries are the same as those observed by Macdonald, Kenderine and Neville 
(1968), who found that the intensity ratios of double- sources in the 3C catalog 
are strongly clustered near unity and that the polarization are strongly 
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clustered around the direction perpendicular to the axis. However, Mitton 
(1972) subsequently found that the latter correlation is less pronounced 
than it had seemed. 

Within this picture, minor deviations from symmetry, which appear in 
many sources, are possible. For example, large scale shearing motions of 
the Intergalactic medium could displace the supercoherent transitions, 
which define the two lobes, in different directions along the field line 
through the central source. This apparent longitudinal displacement of 
the galaxy from the center of the axis joining the lobes could be accompanied 
by a lateral displacement perpendicular to the axis and by unequal emission 
from the lobes. Similar dlstorcions could arise from the uniform motion of 
a rotating galaxy through the intergalactic medium. Another class of 
deviations from the canonical pattern, which obviously fits into the picture 
given here, involves emission from the central galaxy. 

The spatial profile of synchrotron eraissivity, which underlies observed 
maps of radio lr:tensity, is not the same as the profile of electron density 
Fq, for it also depends upon B. More specifically, the emisslvity is given 


by 

c . v-“-« . 


( 87 ) 


where y is the electron spectral index and where the second equality gives 
the dependence expected for a typical radio spectral index of 0.8. In the 
present context, where the theory does not consider perpendicular transport 
and where existing maps do not always resolve structure perpendicular to 
the axis, the most appropriate quantity for comparison with observations is 
the axial profile of power emitted per unit distance parallel to the field. 
From the emisslvity given above and from the relationship « (Fq/B) , it 
follows that this quantity is given by 

(dP/dz) . jlj(v-l) ^ ^0.8 ^-0.8 
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where is the number of particles per unit distance. In the theory 
discussed above, where B depends exponentially upon z, this equation leads 
to a weighting 

(dP/dz) « Fq exp{ -H(y- 1) (z/h) } a: exp{ -0.8 z/L} (89) 

which displaces the radio profile toward the galaxy relative to the density 
profile. Because of this weighting, most of the radio emission comes from 
a tiny fraction of the electrons which are nearest the galaxy in relatively 
strong fields. Moreover, in focused diffusion, the Gaussian peak of the 
radio profile moves toward stronger field'with a velocity given by 

+ (y - 1) (D^^/L) « Y . (90) 

Thus, if a series of discrete explosions occurs in the central galaxy, 
clouds from earlier events drift inward where they 

appear as weak secondary lobes lying on the axis between the intense lobes 
from later events. The following double sources with well marked lobes each 
resolved into a close pair have been reported by Macdonald et. al. (1968) : 
3C33.1, 3C46, 3C61.1, 3C184.1, and 3C234. These authors also report chat, 
in all these cases, the innermost members are weaker than the outermost. 

If the explosions occur frequently or if the acceleration is continuous, 
electrons drifting inward form a continuous bridge between the lobes. 

ClearcuC examples of this behavior are 3C46, 3C274.1, 3C284, 3C430 and 3C452, 
but more or less continuous emission along the axis is seen in many extra- 
galactic sources. Because the Compton/synchrotron mechanism has more time 
to act on the electrons in these inner components^ their radio spectrum is 
expected to be steeper than that of the outer lobes. This prediction is in 
accord with the finding of Macdonald et al. (1968) that the emission from 
between the main components has a steeper spectrum than the source as a whole. 

Radio trails (Wellington, Miley and van der Laan 1973, Miley 1973) are 
believed to delineate magnetic fields dragged out behind as the central galaxy 
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aoves through a stationary medium. In this magnetospherlc configuration 
(Jaffe and Perola 1973), the velocity of the galaxy may exceed the velocity, 
defined by equation (90), with which the radio Gaussian drifts forward. 
Consequently, electron clouds from successive explosions are strung out in 
a trail behind the primary lobes which form, as before, at supercoherent 
transitions moving with the galaxy. This interpretation predicts that 
the radio emission from the tail decreases in intensity systematically with 
distance from the galaxy while its spectral slope increases. These are the 
effects found by Mlley (1973) in 3C129 and NGC 1265. 

Except for the effects interpreted above in terms of ageing and except 
for the flat spectra of central components, which can be similarly inter- 
preted, the structure of radio sources is not strongly dependent upon 
frequency (Macdonald et al.l968, Mackay 1969). This indication that the 
electron spectral index is uniform means that intergalactic propagation is 
not strongly dependent upon rigidity. The same conclusion was reached 
above in regard to interplanetary propagation. 

The time required for electron bunches to reach the supercoherent 
transitions is smaller than that required for the clouds to dissipate. 
Consequently, the probability of observing a double source in its supercoherent 
phase is small. Nevertheless, among the many sources that have been studied, 
it is reasonable to expect that a few are currently in this phase. Such 
sources would appear as two relatively compact radio lobes moving apart with 
a velocity slightly less than twice the speed of light. This superluminous 
velocity of recession is only a little smaller than the velocity found in 
3C279 and is the same as that found in 3C273 by Cohen, et al. (1971, see also 
Whitney, et al. 1971.) Because projection effects lead to an overestimate 
of the velocity of recession when the axis is not perpendicular to the 
observer's line of sight, the predicted velocity seems to be in good sjvree- 
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ment with those observed in these two quasars i 

ConqJton/synchrotron energy loss during the supercoherent phase sets 
an upper limit on the size of radio sources. More specifically, the 
maximum distance z in Kpc at which an electron can arrive with energy E 

11L3X 

in GeV is 

„ -94 X 10^ ^ 2.1 X 10^ (B/v)'" 

[W + .0369 < > ] E W+.0369 <B^> 

-3 2 

where W is the photon energy density in eV cm“^ and < B > is the mean 

square magnetic field in microgauss (pG) averaged over the pitch angle 

distribution and over the distance travelled by the electron at V as c. 

The expression following the second equality, whose form is convenient for 

the analysis of radio data, refers to electrons whose maximum radio 

emission in a field B in pG occurs at frequency v in MHz. If this formula 

is applied to the radio galaxies observed at v = 5000 MHz by Branson et al. 

(1971) and by Pooley and Henbest (1974), who invoked equipartition to deduce 

magnetic fields that can be used to approximate < B > , the calculated values 

of z lie comfortably above the observed semi-major axes by factors of 4 
max 

to 10. Similarly, in the extreme case of 3C236, which was studied at 
V = 612 Iffiz by Willis, Strom and Wilson (1974), where B ~ 0-6 microgauss 
and where the electron lifetime is set by the microwave background (W w 0.38 eV 
cra*^) rather than by the synchrotron effect, z = 16.7 Mpc, which is somewhat 
larger than the distance of 2.8 Mpc between the outer lobes of this giant 
object and the central galaxy. In these examples, energy loss may be 
significant during the evolution of the clouds, but it is insignificant during 
the supercoherent phase. 

In (.ontrast, for the compact outer lobes of Cygnus A (Hargrave and Ryle 
1974), where B » 290 pG and v = 5000 MHz, z^ =16.3 Kpc, which is sub- 
stantially smaller than the 100 Kpc separation of these components 
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from the central galaxy. On the basis of arguments similar to those 

underlying this estimate » these authors reached an equivalent conclusion 

that the synchrotron lifetime of the electrons in these lobes is less than 

the time required for light to reach them for the central galaxy. 

Several considerations soften the impact of this conclusion. To 

obtain more accurate values of z than the crude estimates obtained above 

max 

by substituting in equation (61) the square of the equipartition field, a 

2 

proper evaluation of <B > should take into account not only the collimated 
angular distribution of the supercoherent mode, but also the decrease of B 
with distance from the central galaxy. It is possible that the reduction 
in <B^> arising from the first of these effects significantly outweighs the 
enhancement arising from the second. Of greater significance is the 
possibility that the electron energy density is actually larger than the 
magnetic energy density, for the transient evolution of electron bunches 
and clouds can be controlled by fields weaker than the minimum field 
required for steady state confinement. For the supercoherent phase, this 
possibility seems especially plausible not only because the lateral pressure, 
which tends to disrupt the guiding field, is much smaller for a collimated 
bunch of particles than for an isotropic cloud but also because it is exerted, 
at a given point, for a shorter period of time. Thus, it seems appropriate 
to assume that the compact components are emitted by electrons that have just 
formed superdense clouds after propagating supercoherently across Cygnus A 
in a field of 120 pG corresponding to equipartition within the extended 
components. In this situation, z =61.2 Kpc, which is slightly smaller 

luOX 

than the actual separation of 100 Kpc, but which could probably be made 

2 

consistent by invoking a more accurate value of <B > and by taking into account 
the possibility that equipartition does not apply to the extended components. 
These arguments show that energy loss significantly affects the propagation 
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of electrons in Cygnus A and suggests that the pressure exerted by 
particles may affect the large-scale configuration of the magnetic field. 

If a substantial flux of nuclei accompanies the pure flux of electrons 
assumed above, the latter suggestion becomes conclusive. Thus, outward 
particle pressure, which straightens the lines of force, may underlie 
the well known tendency for the structure of powerful sources like Cygnus A 
to be simpler and more regular than that of weaker sources whose field 
configurations are more easily influenced by motions of the intergalactic 
medium. 

The model of extragalactic radio sources developed here is similar 
to previous interpretations, reviewed by Longair, Ryle and Scheuer (1973), 
which assume that relativistic particles carry energy from the central 
galaxy to the radio lobes. Unlike these interpretations, which invoke 
local acceleration within the lobes, the radio waves are emitted here 
by the same electrons that transport the energy. Fundamental characteristics 
of this transport, which takes place in the same magnetic fields that are 
required to explain the synchrotron emission, give rise to the basic 
morphology of radio sources. Specifically, the emission is confined to 
an axis because electrons propagate parallel to the magnetic field more 
readily than perpendicular to it. Symmetrical lobes appear on this axis 
because electrons are deposited at the supercoherent transitions far from 
the central galaxy where they propagate diffusively. The slow drift 
velocities which characterize this propagation explain the secondary 
structure between the main lobes and establish a relationship between double 
sources and radio trail galaxies. The supercoherent propagation by which 
electrons reach the lobes proceeds at super luminous velocities of recession 
comparable to those observed in some quasars. Thus, focused transport gives 
rise to the radio source structures that are summarized in figure 12. 

Except for relatively slow changes which may occur in powerful sources due 
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to particle pressure, only static magnetic fields are involved. Con- 
sequently, the adiabatic energy losses that embarrass interpretations 
in which electrons are transported within expanding clouds of thermal 
plasma do not occur in the present model. Finally, the basic features 
of focused transport are confirmed by Interplanetary observations. In 
retrospect, it seems surprising that the relationship between the 
scatter— free propagation of solar electrons and the intergalactic 

4 

transport of radio emitting electrons had not been recognized. 
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FIGURE CAPTIONS 

Figure 1. The scattering eigenfunctions at the left, which were 
defined in Paper II, are symmetrical functions of y. In contrast, the 
focusing eigenfunctions at the right are asynouetrlcal. 

Figure 2. The focusing eigenvalues increase monotonically with (V/AL) . 

Figure 3. The four characteristic velocities that reduce in the 
absence of focusing to the finite velocities defined in Paper II all 
vanish when focusing becomes intense. This absence of coupling between 
eigenfunctions of opposite parity leads to supercoherent propagation. 

Figure 4. For isotropic scattering, q = 1, the six characteristic 
velocities that vanish when focusing is absent display a complicated 
dependence upon (V/AL). When the scattering is anisotropic, as It is at 
q = 1.5 and q = 1.9, this intricate pattern is dramatically simplified. 

Figure 5. This graph shows that the artificial decay associated 
with the truncated set of equations Chat describe focused diffusion plays 
an unimportant role in the overall evolution of the distribution function. 

Figure 6. The coherent velocities, and V , that appear in focused 
diffusion are not strongly dependent on (V/AL). 

Figure 7. This graph shows that ~ (1/2L). 

Figure 8. Density profiles which compare ordinary diffusion (a) to 
focused diffusion with positive injection velocity (b) and to focused 
diffusion with negative injection velocity (c). 

Figure 9. The velocities that describe the supercoherent mode are 
plotted as functions of (V/AL). 

Figure 10. A supercoherent density profile is compared to the 
corresponding profile for focused diffusion. 

Figure 11. Schematic diagram of the solar neighborhood showing three 
regions ",n which there appear qualitatively different solar event profiles. 
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Figure 12. An atlas of radio source configurations that can be 
explained in terms of ideas presented in this paper: (a) two symmetrical lobes 

of emission on opposite sides of the central galaxy, (b) two secondary 
lobes between the primary lobes, (c) a succession of lobes deployed 
behind a moving galaxy, and (d) two compact lobes moving away from the central 
galaxy with velocities slightly less than that of light. 
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